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Abstract 

We present a lattice determination of the charm quark's mass, using the mass of the 
Dg meson as experimental input. All errors are under control with the exception of 
the quenched approximation. Setting the scale with Fk = 160 MeV, our final result for 
the renormalization group invariant (RGI) quark mass is Mc = 1.654(45) GeV, which 
translates to m^S(mc) = 1.301(34) GeV for the running mass in the MS scheme. A 6 
percent increase of the RGI quark mass is observed when the scale is set by the nucleon 
mass. This is a typical quenched scale ambiguity, which is reduced to about 3 percent 
for m^^(mc), and to 4 percent for the mass ratio Mc/M^. In contrast, the mass splitting 
rriD* — i^D^ changes from 117(11) MeV to 94(11) MeV, which is significantly smaller than 
the experimental value of 144 MeV. 



1 Introduction 



The charm quark mass is among the fundamental parameters of the Standard Model, and 
its determination from experimental data is of general interest. In practice it is important 
to improve on the current precision: the Particle Data Group in its latest edition gives 
the range 

1.0 < m^(mc) < 1.4 GeV, (1.1) 

thereby increasing the previously quoted uncertainty Q by a factor two. The relatively 
large uncertainty in the charm quark mass has been identified as the dominant error in 
estimates of the fine structure constant at high energies Q, and in phenomenological 
estimates of certain 5-decay rates |QJ5|. 

A determination of quark masses from experiment is necessarily indirect, as quarks 
are not observed as free particles. To establish the connection one needs control over the 
strong interactions at the non-perturbative level. We distinguish two classes of approaches: 
using some variant of the sum rule technique, one relies on perturbation theory and as- 
sumptions such as quark-hadron duality, which allow to connect to quantities derived from 
experiment [P-10|. While the quoted error especially in is rather small, it seems fair to 
say that a reliable assessment of the systematic errors is difficult and may not always be 
possible. 

The lattice approach, on the other hand, is non-perturbative in nature and allows to 
directly compute hadronic observables at fixed lattice spacing and bare quark masses. One 
may then determine the bare quark masses for which the experimental input is matched. In 
order to take the continuum limit one must substitute bare by renormalized masses, which 
are kept fixed as the lattice spacing is varied. It should be emphasized that consistency 
requires the renormalization procedure to be non-perturbative, too. In recent years there 
has been significant progress in non-perturbative renormalization techniques |11]. This has 
led to the determination of the A-parameter and the strange quark mass from low-energy 
hadronic observables, with no uncontrolled systematic errors apart from the quenched 
approximation |12,13]. 

In this paper we present a lattice computation of the charm quark mass, where we 
use the same techniques as for the strange quark mass |13]. In particular we take over the 
results for the non-perturbative quark mass renormalization of ref. |12|, as the renormaliza- 
tion constant is quark mass and hence flavour independent. In addition, all counterterms 
which are needed to remove the leading lattice artefacts in renormalized quark masses have 
been determined non-perturbatively in the relevant range of parameters [|l4,15|. This is an 
important ingredient, as the charm quark mass is not very small compared to the inverse 
lattice spacing, so that cutoff effects can be large. Nevertheless, as will be demonstrated 
below, controlled continuum extrapolations appear to be possible once the leading 0(a) 
artefacts have been cancelled. For this reason our final result for the charm quark mass 
is much more precise than previous lattice estimates [O-p^. In particular, the dominant 



1 



remaining uncertainty stems from the use of the quenched approximation, and further 
progress will require the inclusion of sea quark effects. 

The paper is organized as follows: In section 2 we explain our general strategy and the 
assumptions made. We then describe the technical framework of our calculation and collect 
the formulae needed for the quark mass renormalization and 0(a) improvement (sect. 3). 
Section 4 contains some details of the numerical simulations and a discussion of systematic 
errors. Continuum extrapolations and results for the charm quark mass, the charm to 
strange quark mass ratio and the hyperfine splitting of D mesons are discussed in sect. 5, 
and we conclude in sect. 6. 



2 Strategy 

2.1 Experimental input and basic assumptions 

It is generally assumed that the experimentally observed light hadron spectrum can be 
accounted for by considering the physics of strong and electromagnetic interactions be- 
tween up, down, strange and charm quarks. Furthermore, due to the smallness of the fine 
structure constant, electromagnetic effects are likely to be small. Experimental results for 
the light hadron spectrum can thus directly be compared to results of a QCD calculation, 
possibly after a small correction for electromagnetic effects. 

Numerical simulations of lattice QCD offer the possibility to compute the hadronic 
spectrum with an accuracy of typically a few per cent. One may then turn the tables and 
use the precise experimental data to determine the free parameters of QCD, i.e. the gauge 
coupling and the quark masses. For instance, assuming isospin symmetry, M^^ = M^, the 
four remaining parameters can be determined by matching the kaon's decay constant and 
mass, 

Fk = l60MeV, mx = 495MeV, (2.2) 

the mass of the Dg meson, 

mD, = 1969MeV, (2.3) 
and by taking the quark mass ratio, 

Ms/M = 24.4 ± 1.5, M = i(Mu + Md), (2.4) 

as determined in Chiral Perturbation Theory |21]. This latter input might be traded 



e.g. for the pion mass. However, lattice simulations are typically carried out for quarks 
not very much lighter than the strange quark, so that a chiral extrapolation becomes 
necessary. As the ansatz for such an extrapolation is usually guided by chiral perturbation 
theory, the input of the pion mass is in practice not very different from the direct use of 
eq. (p.4|). A similar problem occurs if one attempts to determine the kaon mass and its 
decay constant at the physical values of the quark masses. Chiral perturbation theory 
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predicts a very weak dependence upon the difference of the valence quark masses [22,2^ 



and this has been verified numerically to some extent |13]. In practice one then computes 
with mass degenerate quarks the sum of which matches the sum of the physical light and 
strange quark masses. Once again, one relies on the assumption that chiral perturbation 
theory provides a reliable description of the quark mass dependence beyond the range 
actually covered by the simulations. 

2.2 Electromagnetic effects 

Electromagnetic effects on hadronic observables involving up, down and strange quarks can 
be computed in Chiral Perturbation Theory. For instance, the result rriK = 495 MeV ( |2.2| ) 
is the (isospin averaged) experimental result after a subtraction of an estimate for the 



electromagnetic mass shift |13|. 



For the charm quark mass determination we would like to estimate electromagnetic 
effects on D-mesons masses, where standard Chiral Perturbation Theory does not applyQ. 
We propose a phenomenological estimate and shall assume that all pseudoscalar D meson 
masses can be parameterized in the form 

mps = AMiight + Bq + C, (2.5) 

where Mught is the light valence quark mass and q distinguishes electrically charged (g = 1) 
and neutral {q = 0) mesons. Note that a linear dependence on the light valence quark 
mass is indeed observed in lattice simulations, at least for some range of light valence 



quark masses [26|. Using this ansatz, we first note that C cancels in the mass differences. 



m^±-m^± = A(Ms-Md), (2.6) 
mj^± - m£,o = A{Ma - Mu) + B. (2.7) 

We may then eliminate A and solve for the electromagnetic mass shift 



B = mD±-moo-j^—j^[mj^±-mo±). (2. 



Taking the mass ratio from Chiral Perturbation Theory |21], 

=(40.3±3.2)-\ (2.9) 

Ms - Md 

and the experimental data for the mass splittings ||2|, 

m^± -mo± = 99.2 ±0.5 MeV, m^± - m£,o = 4.79 ± 0.10 MeV, (2.10) 

^see, however, refs. p4| , p5| ] for an extension of Chiral Perturbation Theory to mesons containing a heavy 
quark. 
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we arrive at the estimate 



B w 2.3 MeV. 



(2.11) 



Similar considerations can be made for the vector mesons, where this effect is even smaher, 
and we conclude that electromagnetic effects cause D-meson mass shifts of a few MeV at 
most. This is below the accuracy currently reached by lattice computations, so that we 
will neglect electromagnetic effects in the following. 

Although there is little doubt that electromagnetic mass shifts are indeed negligible 
we mention that this could also be checked explicitly, by including an additional U(l) 
gauge field in the numerical simulations [27|. 



2.3 Quenched results for the Ajj^^ and Ms 

In the quenched approximation the programme outlined above has already been completed 
except for the charm quark mass, which is the topic of this work. Instead of quoting the 
value of a running coupling in some renormalization scheme at a reference scale, it is 
customary to quote the A parameter in the MS scheme^, viz. [^] 

aJ^ = 238(19) MeV. (2.12) 

Furthermore, the result for the sum of the renormalization group invariant average light 
and strange quark mass is [^] 

M + Ms = 143(5) MeV. (2.13) 

Using the quark mass ratio ( |2.4D , the latter result translates to 

Ms = 138(6) MeV mf^{2 GeV) = 97(4) MeV, (2.14) 

for the running mass in the MS scheme of dimensional regularization at the scale 2 GeV. 
It should be emphasized that these results have been obtained in the continuum limit, 
by carefully disentangling renormalization and cutoff effects. Thus, the only uncontrolled 
error arises from the use of the quenched approximation. 



2.4 The scale tq and the quenched scale ambiguity 

The attempt to describe the real world with data obtained within the quenched approx- 
imation is bound to fail at some point. In particular, the attribution of physical units to 
quenched results is ambiguous. To avoid this ambiguity we will quote results in units of 

^ Note that can be defined beyond perturbation theory, owing to the fact that the exact relation 
between A-parameters of different schemes is determined by the one-loop perturbative relation between 
the respective couplings. Hence, can be defined indirectly, by referring to the A-parameter of a 

non-perturbatively defined renormalization scheme. 
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the low-energy scale rg, which is derived from the force between static color sources [p^ ]. 
This has mostly technical advantages: first, r^ja has been determined very precisely over 
a wide range of cutoff values |2S-^], so that scaling studies are naturally carried out using 



this scale. Second, in the absence of dynamical quarks, tq is only affected by cutoff effects 
of order c? [31|. The relation to other hadronic scales in the quenched approximation is 



known |13|, 



tqFk = 0.415(9), rom^ ^ 2.6. (2.15) 

This illustrates the inconsistency of the quenched approximation, as the quenched result 
for FK/rriN differs by 10 per cent from its experimental value. This can be viewed as a 
scale ambiguity: setting the scale with Fk = 160 MeV is roughly equivalent to the choice 
ro = 0.5 fm, while = 938 MeV corresponds to vq = 0.55 fm. We will later investigate 
the effect of this "quenched scale ambiguity", in order to get a first idea of the size of 
typical quenching errors. 

The results in the previous subsection have been obtained using rg = 0.5 fm. In units 
of ro we have 

r-oAj^ = 0.602(48), roM, = 0.348(13). (2.16) 
2.5 Strategy to compute the charm quark's mass 

At fixed cutoff, standard lattice QCD techniques allow to compute meson masses for given 
bare mass parameters of the meson's valence quarks. As for the strange quark mass, the 



results of [13 1 allow to set its bare mass parameter without further tuning. To determine 
the bare charm quark mass it thus remains to measure pseudoscalar meson masses for 
several bare mass parameters, and to interpolate the meson masses versus the bare charm 
quark mass to the physical point where 

mz)^ = 1969 MeV romu, = 4.99. (2.17) 

Then, using known renormalization factors, one may obtain the renormalization group 
invariant charm quark mass, Mc, in units of tq and at the given value of the cutoff o~^. 
Repeating the procedure for smaller lattice spacings a, one may eventually extrapolate to 
the continuum. Additional control of this extrapolation may be obtained by considering 
alternative definitions of vqMc, with different cutoff effects but the same continuum limit. 

In the following two sections we describe in some detail the technical setup to achieve 
this goal. The reader not so much interested in the technicalities may directly go to sect. 5, 
where the continuum extrapolations and final results are discussed. 

3 The technical framework 

For the practical implementation of the programme we use the framework of 0(a) improved 



lattice QCD as described in ref.|32|. For unexplained notation we refer to this paper. 
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3.1 SF correlation functions 



In order to extract hadron masses we use correlation functions derived from the QCD 



Schrodinger functional (SF) |^,34|. Hence, QCD is considered on a space-time manifold 
which is a hyper-cylinder of volume x T, with the quantum fields satisfying periodic 
boundary conditions in space, and Dirichlet boundary conditions at Euclidean times xq = 
and xq = T. The technical advantages for the computation of hadron properties have 
been demonstrated in |^^. One of the main points is the possibility to define quark and 
antiquark boundary states C and which only transform under spatially constant gauge 
transformations. Boundary sources with quantum numbers of mesons are therefore gauge 
invariant even if quark and antiquark fields are localized at different points in space. 

For the study of pseudoscalar and vector mesons, we define the axial current and 
density, 

Af,{x) = i)i{x)j^j5ll)j{x), P{x) ='ipi{x)j5lljj{x), (3.18) 

and the local vector current, 

V^{x) =i^i{x)j^ij,{x), (3.19) 
with flavour indices i ^ j- The correlation functions 

fA{xo) = -^a6^(C,(y)75Ci(z)Ao(x)>, (3.20) 

/p(^o) = -^a6^(0(y)75Ci(z)P(x)>, (3.21) 

kyixo) = -^a6^(0(y)7kCi(z)Vk(x)), (3.22) 
y>z 

are then used to compute pseudoscalar and vector meson masses with valence quark 
flavours i and j. The notation and technology has been described in detail in ref. [p^ ]. 
One studies the effective masses 

{xo + ia) = In {/(xo)//(xo + a)} , (3.23) 

where /(xq) stands for one of the above correlation functions. For large time distances 
from the boundaries, i.e. large xq and T — xq, one expects the ground state of a given 
channel to dominate the correlation function, which manifests itself in a plateau for the 
effective mass versus xq. 

3.2 Definition of renormalized quark masses 

The renormalization of quark masses for Wilson type quarks is complicated by the fact 
that all axial symmetries are explicitly broken by the regularization. As a consequence. 
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quark mass renormalization is both additive and multiplicative, the axial current requires 
a scale-independent renormalization, and axial Ward identities, such as the PCAC relation 
are violated by cutoff effects. There are various ways to define renormalized quark masses, 
which are equivalent in the continuum limit, but may differ at finite lattice spacing. This 
will later be used to to achieve a better control of the continuum extrapolations. 

We start with the definition of a bare current quark mass using the PCAC relation, 



2/p(xo) ' 



m,^ = ^^^l^, (3.24) 



where 9^ = ^{d^ + d*) denotes the symmetric lattice derivative in ^u-direction. Prom the 
bare quark mass niij we obtain the sum of the renormalized valence quark masses through 
multiplicative renormalization of the axial current and density in ( |3.24D , viz. 

rriR^i + mjij = 2ZxZ^^mij. (3.25) 

A single quark mass is obtained if the quark flavours i and j are mass degenerate, 

W'R.i = Zp^Z^^rrii. (3.26) 

Another starting point is provided by the hopping parameter in the Wilson quark 
action. One defines a bare subtracted quark mass, 

"lq,i = ^(N"^ - '^c~ritical)' (3-27) 

where ^critical is the value of k at which chiral symmetry is restored. This subtracted bare 
quark mass is then multiplicatively renormalized, 

m.R,i = .Zm"lq,i- (3.28) 

In what follows, we will always refer to the same renormalization scheme for the quark 
mass. This is achieved by first defining the renormalized axial current and density, and 
by setting 

Z^ = ZZ^/Zp, (3.29) 
where Z is the scale independent ratio between the bare quark masses, 

nii = Zniq^i. (3.30) 

Furthermore, we restrict attention to quark mass independent renormalization schemes. 
The SF scheme of ref. p6| has this property, and it has the additional virtue of being 
defined beyond perturbation theory. Its relation to the renormalization group invariant 
(RGI) quark mass, 

M, = lim [2bogHfi)] , (3.31) 
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has been determined for A'^f = ||T2|. We recall that bo = (11 - |A'f)/(47r)^ and do = 
8/(47r)^ are the lowest order perturbative coefficients of the renormalization group func- 
tions for gauge group SU(3). The running mass coincides with at some scale 
fj, = fj,Q and is otherwise determined by the RG equations. The RGI quark mass is non- 
perturbatively defined, it is scheme and scale independent and hence a natural candidate 
for a fundamental parameter of QCD. However, it is still customary to quote quark masses 
in the MS scheme of dimensional regularization at some reference scale. Once M is given, 
the ratio m^^{n)/M can be computed by using the perturbative renormalization group 
functions, which are known with 4- loop accuracy [07|j3^]. 



3.3 0(a) improvement 

In practice, cutoff effects with Wilson quarks can be rather large, and it is advisable to 
cancel at least the leading 0(a) artefacts in physical observables. One may distinguish 
cutoff effects which arise already in the chiral limit and those proportional to the quark 
masses. To cancel the former, a single counterterm, the so-called Sheikholeslami-Wohlert 
term ||39[] , must be included in the action, with a coefficient Csw which is a function of the 
bare coupling only. For the improvement of on-shell correlation function one also needs 
to improve the operators. In our context, only the improvement of the axial current is 
relevant, 

(^i)^ = A^ + CAad^P, (3.32) 

with the improvement coefficient ca. Both coefficients are known non-perturbatively for 
A^'f = and bare couphngs 5o = < 1 ||40t| . 

0(a) artefacts which are proportional to the quark mass ruq can be dealt with per- 
turbatively, as long as aniq is small enough. In practice this may still be the case for the 
strange quark, but in the case of charm the bare quark mass in lattice units is typically in 
the range 0.25 — 0.5. In the quenched approximation, the improvement coefficients needed 
to compute 0(a) improved renormalized quark masses are all known non-perturbatively. 
Renormalization and 0(a) improvement of the subtracted bare quark mass is achieved 
separately for each quark flavour through 

m-R = Zmmq{l + bmamq). (3.33) 

As for the PCAC quark mass, we ffi'st recall the form of the renormalized 0(a) improved 
axial density and current, 

(Ar)^ = Za [l + 6A^(amq,j amqj)] (yli)^, (3.34) 
Pr = Zp [l + bp^{amq^i + aniqj)] P. (3.35) 

For the quark mass improvement this means 

w-R,i + "^RJ = 2ZaZp^ [l + (6a - bp)^{amq^i + arnqj)] nnj. (3.36) 
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Here, it is implicitly assumed that the bare mass rriij is defined including the counterterm 
cx CA, i.e. 

_ dpfAjxp) + CAa9o9o/p(xo) 

Besides this definition with standard lattice derivatives we also use improved derivatives 
as in ref. While the difference is formally an effect of O(a^) it may be important 

numerically. 

The improvement coefficients and &a — bp have been determined for A'^f = in the 



relevant range of bare couplings [14|. This, together with the renormalization constants 
Z 1 14], Za [41 1 and Zp [^] allows us to define the renormalized 0(a) improved charm 
quark mass in the SF scheme in various ways. One may then convert from the SF scheme 
to the RGI quark mass, using the flavour independent ratio M/m(//o) in the continuum 



limit 1 12]. It is convenient to combine the factors relating the bare current quark mass 



and the RGI masses into a single one, 

M 



m(/xo) 



ZA^pi, (3.38) 



which has been parameterized in ]0] over some range of bare couplings. Up to O(a^) we 
then have 

Mi = Zm [1 + (^A - bp)am^.i)\ mi = ZMZm^.i{l + bj^amq^i), (3.39) 
and an analogous equation for Afj + Mj in the mass non-degenerate case. 

4 Numerical simulations 

The numerical simulations were carried out on machines of the APEIOO and APEIOOO 



series using single precision arithmetic and standard algorithms (see ref. ]40| for a detailed 
description and references). Autocorrelation times were found to be small so that our 
measurements could be treated as statistically independent for all observables considered. 

4.1 Simulation parameters 

Our choice of simulation parameters is displayed in table |l[ The non-perturbatively de- 
termined coefficients Cgw and ca were taken from ref. 1^0|. We chose the same four values 
of P as in [O], which correspond lattice spacings in the range a = 0.05 — 0.1 fm. The con- 



version to physical units is done using the scale rg = 0.5 fm j28] and the parameterization 



of ref. H, 



ln(a/ro) = -1.6805 - 1.7139(/3 - 6) + 0.8155(/3 - 6)^ - 0.6667(/3 - 6)^ (4.40) 
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The lattice size Lja varied between 16 and 32, such that the hnear extent of the spatial 
volume was around L ~ 1.5 fm or larger. For spatial volumes of this size we do not expect 
any sizeable finite volume effects on our pseudoscalar and vector masses. This is supported 
by previous findings e.g. in ref. where the pseudoscalar meson mass was found to be 
safe against finite volume effects, provided its Compton wave length was about 5 times 
smaller than L. Similar statements are expected to hold for the vector mesons. In the 
time direction the lattice size varied between 40 and 80, such that the physical extent T/rg 
was kept constant to a good approximation. 



/? 


Tja 


Lja 




T/ro 


ro/a 


-^mcas 


6.0 


40 


16 


2.98(1) 


7.45(3) 


5.368(22) 


350 


6.1 


48 


24 


3.80(2) 


7.59(3) 


6.324(28) 


94 


6.2 


54 


24 


3.26(2) 


7.34(4) 


7.360(35) 


159 


6.45 


80 


32 


3.06(2) 


7.65(4) 


10.458(58) 


123 



Table 1: Simulation parameters, lattice sizes and number of independent measurements. 
Our estimates of the hopping parameter for the strange quark are based on table 2 



of ref. 1 13 1 where the quantity 

[ZMi?/ro],2„2^=i,5736 = 1.5736 x r^{M, + M), (4.41) 

is quoted for each value of /?. In order to isolate the strange quark mass we use the mass 
ratio from chiral perturbation theory ( |2.4| ). Then we relate the RGI strange quark mass 
to the subtracted bare quark mass in lattice units, 

roMs = (ro/a)ZAfZamq,s(l + 6mamq,s). (4.42) 



where we use the coefficients Z and 6m of ref. jlj] . Finally, to relate to the corresponding 
hopping parameter Kg ( 3.27] ), we use estimates of ^critical at /3 = 6.0,6.2 and inter- 



polations at /3 = 6.1,6.45 (cf. table |^. In this way the hopping parameter for the 
strange quark could be fixed prior to the simulations. As usual, there is some uncertainty 
associated with the choice of ^critical- In order to assess the sensitivity of our observables 
to the precise choice of the strange quark mass, we also used a second value of Kg for all 
but our largest /3- value. Finally, we chose 3 k- values around the expected physical charm 
quark mass. The numerical values for all simulated hopping parameters are displayed in 
table 0. 



4.2 £)-meson masses 

From the correlation functions ( 3.20| - O^ ) we determined the corresponding effective 
masses myjyXi^ in the pseudo-scalar (X = PS) and vector channel (X = V). In addition we 
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K 


aruq 


am 




6.0 


0.1190 
0.1200 
0.1210 
134108 
0.133929 


0.5033 
0.4683 
0.4339 
0300 
0.0350 


0.5129(2) 
0.4736(2) 
0.4357(2) 
0304f3) 
0.0353(3) 


0.2962(3) 
0.2913(3) 
0.2837(2) 
0299f4) 
0.0346(4) 


6.1 


0.1218 
0.1224 
0.1230 
134548 
0.134439 


0.4149 
0.3948 
0.3749 
0260 
0.0290 


0.4425(2) 
0.4197(2) 
0.3974(2) 
0317(3) 
0.0348(3) 


0.3113(2) 
0.3030(2) 
0.2940(2) 
0312f3) 
0.0342(3) 


6.2 


0.1230 
0.1235 

U.iz4U 

0.134959 
0.134832 


0.3830 
0.3666 

U.ooUz 

0.0228 
0.0263 


0.4072(2) 
0.3888(2) 
U.o7UD(zj 
0.0254(2) 
0.0291(2) 


0.3086(2) 
0.3000(2) 
0.2909(2 j 

0.0253(3) 
0.0288(2) 


6.45 


0.1270 
0.1280 

0.1290 
0.135124 


0.2524 
0.2217 
0.1914 
0.0157 


0.2646(1) 
0.2323(1) 

0.2006(1) 
0.0166(1) 


0.2336(1) 
0.2094(1) 

0.1844(1) 
0.0165(1) 



Table 2: Simulated K-values and corresponding results for the bare subtracted and bare 
current quark masses from the mass degenerate PCAC relation. 

considered the effective mass ms associated with the ratio /{xq) = kylxQ)/ fp{xo), which 
directly yields the mass splitting between the vector and the pseudoscalar channel. For all 
parameter choices these effective masses exhibit a plateau. Deviations from the plateau 
value at small xq or small T — xq are expected, due to the contribution of excited states. 
Furthermore, since the masses in our work are quite large, the correlation functions decay 
rapidly and one might expect problems with rounding errors at the larger values of xq. 

4.2.1 Rounding errors 

To estimate rounding errors, we varied the precision of the quark propagator compu- 
tation. In our production runs the Wilson-Dirac equation was solved with a squared 
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Analysis of rounding errors 




2 4 6 8 

Figure 1: Illustration of the rounding error problem: The plot shows, for all 4 /?- values, the 
logarithm of the relative differences between two effective masses which differ by the solver 
precision used (cf. text). The curves are obtained as averages over 0(10) configurations 
using the Kc-value closest to the physical point. The chosen plateau region is indicated by 
the vertical lines. 

relative precision = 10"^^. In addition, at each /? value we chose ten independent 
gauge field configurations, and repeated the calculation with = 10~^^ and = 10~^^. 
Indicating the solver precision by a superscript, we find that the ratio of effective masses 
myS'^'^^ /myS'^'^^ — 1 grows roughly exponentially as a function of the time coordinate. We 
performed an exponential fit to this ratio of the form rexp(xo-R), which was then used to 
estimate rounding errors on the effective masses. Requiring that the relative error on the 
masses not exceed a certain value (typically a fraction of a percent) then leads to an upper 
limit for the range of the plateau. However, note that our method is likely to overestimate 
this systematic effect, as the limit of single precision is only reached for solver residuals 
between = lO'^^ ^nd = 10"^^. 
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Figure 2: The effective mass plateaux from the correlation function fA{xo) for 3 charm 
quark mass parameters at all /3 values. The plateaux regions are indicated by the vertical 
lines. The pseudoscalar meson masses cover the physical Dg meson mass (full line) except 
at 13 = 6.2. 

4.2.2 Excited states 

To estimate the contributions of excited states to the effective masses we first chose an 
a priori range for the plateau. Then we averaged nixixQ) over the plateau region and 
subtract the averaged value from mx(a^o)- For small xq we fit the remainder with an ansatz 
of the form r/exp(— xqA). Reasonable fits could be obtained at all (3 values, and could 
be used to quantify the contamination by excited states from the lower time boundary. 
Prom the upper boundary at xq = T one expects a contribution from the O'*'^ glueball 
with mass roiriQ ~ 4.3 [^5|. While in some cases a rough confirmation seemed possible, we 
essentially assumed this to be the case. In fact, given the problem with rounding errors 
at large values of xq, a clean signal for the glueball state would have come as a surprise. 
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4.2.3 Plateau regions 



Adding all the systematic errors linearly, a plateau region was determined as the interval 
of times xq where the relative systematic error was smaller than some prescribed value. 
We chose 0.3% for the pseudoscalar mass, 0.5% for the vector mass and 0.8% for the mass 
splitting. The plateau ranges found in this way were remarkably stable when expressed 
in physical units. Therefore we decided to define the final ranges as follows: from 3.8 
to 5.0 ro for the pseudoscalar masses, from 3.8 tq to 5.2 rg for the vector masses and from 
3.7 ro to 4.5 ro in the case of the mass splitting. At fixed f3 and for fixed quark mass 
parameters, the meson masses are then obtained as averages over the plateau region. 
Statistical errors were determined by a jackknife procedure, and the maximally allowed 
systematic error was taken as final systematic error on the effective masses. This is likely 
to be an overestimate, as the systematic error on the effective masses inside the plateaux 
is a bit smaller. The results obtained in this way are displayed in table |5|. 



4.3 Determination of the charm quark mass 

The meson masses at the simulated parameter values, can be considered as functions of 
the corresponding bare valence quark masses. Besides the bare subtracted quark masses 
(obtained using the ^critical values in table we used bare PC AC masses from the SF cor- 
relation functions at xq = T/2. An analysis of rounding errors for the PCAC masses shows 
that these are at the level of 0.01% and hence completely negligible. We distinguish the 
PCAC relation with degenerate and non-degenerate quark masses. Furthermore, we con- 
sider two choices of the lattice derivatives used in the PCAC relations: first, the standard 
choice involving differences between nearest neighbours only, and, secondly, the improved 
derivatives which also involve next-to-nearest neighbours [cf. [|l^] for the definitions]. The 
results are given in tables |2| and |5|. 

We recall that the strange quark mass parameters were already chosen according 
to eq. (4.42). Moreover, it turns out that the meson masses only mildly depend on the 



strange quark mass. Hence it remains to determine the bare charm quark masses where 
the pseudoscalar meson mass assumes its experimental value, 

romps = roniDs = 0.5 fm x 1969 MeV = 4.99. (4.43) 

Our simulation parameters are such that this point could be reached by an interpolation, 
with the exception of f3 = 6.2, where a small extrapolation was required. We performed 
linear fits of the form 

rompg = ao + airom, (4.44) 

where m is some definition of the bare charm quark mass. The fits were done inside a 
Jackknife procedure to take into account correlations of the data points. 
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Figure 3: The inter- or extrapolation of mps as determined from figure 2 vs. the bare 
charm quark mass rorrisc- 

Given the various interpolated bare charm quark masses we construct the correspond- 
ing 0(a) improved RGI masses. We distinguish 5 definitions, the first 3 being given by 

'^o^cU.c = ^A/|2romsc [l + (6a - 6p)^(a"iq,c + aruq^s)] 
- rorris [1 + (6a - 6p)amq,s)]|, 

roMclm, = ZMrorricil + {bA - bp)aniq^c] , 
roMclm^^^ = ZM^'romq^c [1 + 6mamq,c] • (4.45) 

In addition we use the analogues of the first two definitions but with mgc — > tusc"^ and 
rric — > rn^^, i.e. using the next-to-nearest derivatives in the PCAC relations. The results 
are shown in table ^, where the total error is obtained by including the errors of ro/a, and 
of the renormalization constants and improvement coefficients in quadrature [cf. table 
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^ '^critical 'Z'M 



Z bA-bp bra 



6.0 0.135196 1.752(19) 1.0604(4) 0.171(5) -0.709(6) 

6.1 0.135496 1.782(20) 1.0852(5) 0.071(3) -0.699(5) 

6.2 0.135795 1.807(20) 1.0960(5) 0.039(3) -0.691(7) 
6.45 0.135701 1.852(20) 1.1045(5) 0.010(5) -0.673(12) 

Table 3: Renormalization constants and improvement coefficients at the simulated (5- 
values. 



/3 




roMc\m^ 


?'oM;|mq,c 




roMc\m)}^ 


romy 




6.0 
6.1 
6.2 
6.45 


4.331(59) 
4.274(59) 
4.277(55) 
4.220(60) 


5.215(75) 
4.824(70) 
4.682(67) 
4.428(64) 


3.224(41) 
3.479(43) 
3.711(47) 
3.975(53) 


3.566(47) 
3.755(50) 
3.905(52) 
4.038(57) 


3.026(37) 
3.492(44) 
3.769(49) 
3.995(56) 


5.281(14) 
5.299(15) 
5.299(15) 
5.300(22) 


0.299(11) 
0.303(12) 
0.299(13) 
0.300(17) 


CL 


4.19(11) 


4.20(12) 


4.27(10) 


4.21(11) 


4.31(11) 


5.300(35) 


0.297(26) 



Table 4: Results for the five definitions of the RGI charm quark mass [cf. eqs. ( 4.45| )], the 
vector meson mass and the mass splitting between vector and pseudoscalar mesons. 

5 Continuum extrapolations and results 
5.1 RGI charm quark masses 

We now come to our main results, the continuum extrapolation of the RGI charm quark 
mass. As discussed in sect. 3, the charm quark mass can be defined in various ways 
which are all equivalent up to cutoff effects. Since complete 0(a) improvement has been 
implemented, we attempt linear fits to the data of the form 



rQM,=A + B{a^/rl), 



(5.46) 



with the two parameters A and B. As can be seen in fig. 4, the fits to the data appear very 
reasonable. Excluding the coarsest lattice spacing, the per degree of freedom ranges 
from 0.1 to 1.5. Moreover, the various definitions of vqMc all yield compatible continuum 
extrapolated results. We also performed combined fits, with the constraint of a common 
continuum limit. The correlation of the data at the same /3-value was taken into account 
by determining the covariance matrix in a Jackknife procedure. However, the correlation is 
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Figure 4: The continuum extrapolation for all 5 definitions of vqMc considered. In the case 
of roMcl^i^m higher order cutoff effects become visible. For illustration we also included a 
fit to all 4 data points allowing for an additional term oc (a/ro)^ (dashed line). 



rather strong, so that the error is not much reduced compared to the individual continuum 
extrapolations. 

As our best result we choose the one obtained from the heavy-light PCAC relation, 
as cutoff effects are found to be rather small in this case. We quote 

roM, = 4.19(11), (5.47) 

where the error contains the 1.3 percent error of the matching factor to the RGI mass []l^ ]. 
In physical units we then have (with tq = 0.5 fm) 

Mc = 1.654(45) GeV. (5.48) 



5.2 Conversion to the MS scheme 

In our framework the RGI quark masses appear as the primary quantities. However, many 
phenomenological applications use the masses in the MS scheme at the scale set by the 
quark mass itself. In order to convert to this scheme, we integrate the renormalization 
group equations using the known 4-loop KG functions for A'^f = |37,^|, and find 



mf ^(mc) = 1.301(28)(20) GeV = 1.301(34) GeV, (5.49) 

where we have used the result ( ^.161) for A^. The two errors given in the first number 
correspond to the error in the mass and the change induced by a variation of a|^ 
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within its error bars respectively. As both errors are independent we have combined them 
quadratically in the last equation of ( 5.4S| ). Note that the scale of the charm quark mass 



is already quite low, so that the order of perturbation theory must be specified, too. In 
fact, using 3 and 2-loop RG evolution we obtain, 

_MS,-^ /l.294(34)GeV 3-loop evolution, 

(mc) = < (5.50) 
[1.263(34) GeV 2-loop evolution, 

i.e. even the 4-loop contribution is still sizeable, and the difference between 2-loop and 
4-loop RG evolution is as large as the total error. 



5.3 Further results 

5.3.1 The mass ratio Mc/M^ 

A further consistency check for our results is provided by considering the ratio Mc/Mg. 
Here one may compare the continuum extrapolation of the ratio Mc/Mg, to the ratio of the 
two continuum extrapolated quark masses. Unfortunately, this exercise requires a rather 
precise tuning of the strange quark mass. As our strange quark mass parameters were 
tuned using interpolated values of ^critical) the corresponding curve Ks(/3) does not very 
precisely correspond to a condition of constant physics. Therefore, instead of using our 
own data for the strange quark mass, we decided to use the bare current quark masses from 



ref. [13,^], where a careful extrapolation to the physical kaon mass has been performed. 
After 0(a) improvement of the bare masses, we take the ratios and combine the errors in 
quadrature. The continuum extrapolation then yields 

Mc/Ms = 12.0(5), (5.51) 

which can be compared to the ratio taken directly in the continuum limit, viz. 

Mc/Ms = 12.2(1.0). (5.52) 



5.3.2 n^D* and rriDs ~ ''^D* 

We also computed the mass of the vector meson D*. The cutoff effects are very small, so 
that the continuum extrapolation is not problematic. We obtain 

roniD* = 5.300(35) ^ mn* = 2092(14)MeV. (5.53) 

Despite the quenched approximation this is not far from the experimental result, mu* = 
2112 MeV ||2[. It is customary to study the mass splitting ms = tud* — m£)^. Compared to 
the experimental value of 144 MeV the results of quenched lattice simulations often turn 
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out to be smaller (see, e.g. |26|). The mass splitting may be obtained in two ways: we 
may either subtract the input value romog = 4.99 from the result ( p.53| ), which yields 



roruD- - 4.99 = 0.310(35) = 122(14) MeV x 0.5 fm. (5.54) 

On the other hand, the mass splitting may be directly obtained by studying the effective 
mass associated to the ratio of correlators /(xq) = ky{xQ)/ fp{xo) (cf. sect. 2). With this 
direct method we obtain in the continuum limit 

roiniD* -moj = 0.297(26) = 117(11) MeV x 0.5 fm, (5.55) 

which agrees with the previous result within errors. Both continuum extrapolations are 
shown in figure 5, and we conclude that the mass splitting is indeed smaller than the 
experimental value, albeit only by 1.5 and 2.5 standard deviations for the indirect and 
direct methods respectively. 



5.4 Quenched scale ambiguities 

As mentioned in sect. 2, the quenched approximation to QCD accounts for the observed 
hadronic spectrum up to inconsistencies which are at the 10 per cent level. Turning this 
around, the choice of different hadronic input from experiment will lead to a spread of 
results for the A parameter and the quark masses in the quenched approximation. Under 
the assumption that these inconsistencies are entirely due to the neglected sea quark 
effects, this spread of results may be taken as a first estimate of the quenching error. 

We start with some definition of the RGI charm quark mass, consider it as a function 
oi z = roniDs , 

roM, = F{z). (5.56) 
and expand around the standard choice used so far, zq = 4.99, 

F{z) = F{zo) + {z- zo)F'{zo) + O {{z - zof) • (5.57) 

For a 10% shift in the scale rg we have z — zo = 0.5. With this choice we obtain an estimate 
of the first order term in the continuum limit, 

0.5F'(4.99) = 0.7(1). (5.58) 

Expressing Mc again in physical units (using now ro = 0.55fm!) we observe that the 
charm quark mass increases by about 6 percent. Although this number is only considered 
a rough estimate we remark that the approximations made are supported by our data: 



first, figure 3 indicates that the higher order terms in eq. (5.57) are indeed small over quite 



some range of charm quark masses. Second, we have assumed that the shift in the strange 
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Figure 5: Results for the mass splitting mg = m^* —rrij:)^. The open points are obtained by 
computing the vector meson mass while the solid points represent the direct computation 
of mg from a ratio of correlation functions (the data points have been slightly offset for 
clarity) . 



quark mass can be neglected. As observed in the relative change in the light quark 
masses is the same as for ro itself, due to the proportionality 



M oc rompg. 



(5.59) 



As already mentioned earlier the sensitivity of the Dg meson mass to a variation of the 
strange quark mass is rather low, and this can be seen in table ^, where the differences 
between the two chosen Kg-values correspond to 10 — 15 percent differences in the quark 
mass. 

Having established the 6 percent increase of Mc, we note that this corresponds to 



an increase of only 3 percent for m, 



This is due to the quark mass anomalous 



dimension in the MS scheme, which is such that an increase of the renormalization scale 
decreases the running quark mass. In passing we mention that this also holds for the 
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strange quark mass at the scale fi = 2 GeV where the induced change amounts to about 8 
percent. 

The same exercise can be done for the vector meson mass and the mass sphtting. The 
first order terms are now given by 

f 0.46(1) foT F = roruD* , 
0.5F'(4.99) = < " (5.60) 

1—0.035(5) for F = ro(m£)* — m£)^). 

When expressed in physical units this corresponds to changes of +9% and —20% respec- 
tively. Hence, the quenched scale ambiguity for the mass splitting is rather large. In 
particular, in this latter case the discrepancy to the experimental result increases to about 
4 standard deviations. 

Finally we look at the ratio Mc/Ms, which may be expected to be less sensitive to 
the scale ambiguity. Here, the 6% shift of the charm quark mass is overcompensated by 
the 10% shift in the strange quark mass, leaving a quenched scale ambiguity of about 4%. 



5.5 Taking m/j. as input 

Finally we consider setting the charm quark mass using the L'*-meson instead of Ds 
as input. The analysis is completely analogous and the continuum extrapolations look 
qualitatively very similar. The Dg mass becomes now a measurement, and the result, 

roDs = 5.04(4), (5.61) 

is not far from the experimental value 4.99 previously used as input. As our best result for 
the charm quark mass we quote again the value obtained from the mass non-degenerate 
PC AC mass. The result, 

roMc = 4.27(13), (5.62) 



is slightly larger but consistent within errors with eq. (5.47). 



6 Conclusions 

The main result of this paper is the determination of the RGI charm quark mass using 
the experimentally measured Ds meson mass as essential input. Apart from the quenched 
approximation all errors appear well under control. With specified input the total error 
is around 3 percent, which is smaller than the error one might attribute to the use of 
the quenched approximation. An estimate of the latter has been obtained by varying the 
scale which is used to assign physical units; assuming that the real world is described by 
full QCD, the inconsistencies may be taken as a first indication of the quenching error. 
Surprisingly, a scale variation of 10% induces only a 6% shift of the RGI charm quark 
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mass, which is further reduced to 3% for the charm quark mass in the MS scheme at the 
charm quark mass scale. 

While the agreement with some non-lattice charm quark mass determinations is sur- 
prisingly good (in particular the value m^^(mc) = 1304(27) MeV of ref. 0) we emphasize 
that the real quenching error can only be asserted by going beyond the quenched approx- 
imation. However, our results show the way towards more realistic studies, as well as the 
potential strength of lattice techniques. In particular, we draw the following conclusions: 

• precise results in charm physics are attainable using the standard set-up of 0(a) 
improved lattice QCD, which was originally designed for the light quarks, 

• cutoff effects can be quite large in general, making a continuum extrapolation neces- 
sary. In our examples decent continuum extrapolations were possible based on data 
covering a factor of 2 in the cutoff scale. 

Finally, as a by-product of our simulations we measured the D* meson mass and the mass 
splitting between the vector and pseudoscalar Dg mesons. The latter turns out to be 
rather small when compared to experiment, and is very sensitive to the quenched scale 
ambiguity. It will be interesting to look at the effect of sea quarks on the mass splitting. 

This work is part of the ALPHA collaboration research programme, and partially sup- 
ported by the European Community under the grant HPRN-CT-2000-00145 Hadrons/Lattice 
QCD. Simulations were carried out on machines of the APEIOO and APEIOOO series at 
DESY-Zeuthen. We thank the staff at the computer centre for their help, and P. Ball, 
B. Bunk, T. Hurth, R. Sommer, H. Wittig and U. Wolff for useful discussions. 
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A Raw lattice data 









arusc 


imp 

anise 


ampg 


amy 


arrig x 10^ 


6.0 


0.1190 


0.134108 


0.2356(3) 


0.1962(4) 


0.9306(16) 


0.985(3) 


5.55(20) 






0.133929 


0.2387(3) 


0.1982(4) 


0.9366(15) 


0.990(3) 


5.43(20) 




0.1200 


0.134108 


0.2214(3) 


0.1869(4) 


0.8998(16) 


0.957(3) 


5.80(20) 
















O.O t yzyj J 




0.1210 


0.134108 


0.2074(3) 


0.1775(3) 


0.8683(15) 


0.928(3) 


6.06(21) 






0.133929 


0.2104(3) 


0.1796(4) 


0.8744(15) 


0.933(3) 


5.93(20) 


6.1 


0.1218 


0.134548 


0.2125(3) 


0.1871(3) 


0.811(2) 


0.858(3) 


4.68(21) 






0.134439 


0.2144(3) 


0.1885(3) 


0.815(2) 


0.862(3) 


4.64(20) 




0.1224 


0.134548 


0.2038(3) 


0.1806(3) 


0.792(2) 


0.841(3) 


4.81(21) 












u. / J 








0.1230 


0.134548 


0.1951(3) 


0.1740(3) 


0.772(2) 


0.823(3) 


4.96(21) 






0.134439 


0.1970(3) 


0.1754(3) 


0.776(2) 


0.826(3) 


4.92(20) 


6.2 


0.1230 


0.134959 


0.1958(2) 


0.1770(2) 


0.730(2) 


0.769(3) 


3.73(19) 






0.134832 


0.1980(2) 


0.1787(3) 


0.735(1) 


0.773(2) 


3.67(17) 




0.1235 


0.134959 


0.1886(2) 


0.1713(2) 


0.714(2) 


0.754(3) 


3.84(18) 






0.134832 


0.1907(2) 


0.1730(2) 


0.719(1) 


0.758(2) 


3.77(17) 




0.1240 


0.134959 


0.1814(2) 


0.1655(3) 


0.698(1) 


0.739(3) 


3.95(19) 






0.134832 


0.1835(2) 


0.1673(3) 


0.703(1) 


0.743(2) 


3.88(17) 


6.45 


0.1270 


0.135124 


0.1332(2) 


0.1270(2) 


0.518(2) 


0.544(3) 


2.59(18) 




0.1280 




0.1191(2) 


0.1142(2) 


0.482(1) 


0.512(3) 


2.82(20) 




0.1290 




0.1049(2) 


0.1012(2) 


0.446(1) 


0.478(3) 


3.09(20) 



Table 5: Results for unrenormalized current quark and meson masses for all heavy-light 
combinations of simulated k- values (errors are statistical only). 
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